Discontinuous Absorbing State Transition in (1 + 1) Dimension 
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A (1 + 1) dimensional model of directed percolation is introduced where sites on a tilted square 
lattice are connected to their neighbours by N channels, operated at both ends by valves which are 
either open or closed. The spreading fluid is assumed to propagate from any site to the neighbours 
in a specified direction only through those channels which have open valves at both sites. We show 
that the system undergoes a discontinuous absorbing state transition in the large N limit when the 
number of open valves at each site n crosses a threshold value n c — y~N. Remarkable dynamical 
properties of discontinuous transitions, like hysteresis and existence of two well separated fluctuating 
phases near the critical point are also observed. The transition is found to be discontinuous in all 
(d + 1) dimensions. 
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Models of directed percolation (DP)[1] mimics filtering 
of fluids through porous media along a given direction. 
These models show a continuous phase transition from 
a non-percolating (absorbing) to a percolating (active) 
state when connectivity of pores is increased. Numer- 
ous physical phenomena, like forest firedJ], epidemics Q, 
transport in random mediaQ, synchronization^ and 
chromatography can be understood using these models. 
The corresponding critical behaviour forms a universality 
class of absorbing state phase transitions, formally known 
as the DP-class, which has been realized [6] convincingly 
in (2+1) dimension. 

In (1 + 1) dimension, the directed bond percolation[7( is 
modelled as a spreading process on a tilted square lattice 
where bonds connect any two neighbouring sites with 
probability p. A bond represents an open channel con- 
necting two pores. On this random structure, the spread- 
ing agent (or the fluid) percolates along a given direction 
through the bonds (open channels). The system is said 
to be percolating if the agent spreads across an infinite 
cluster of connected bonds with finite probability. It is 
well known that percolation occurs here as a continuous 
(non-equilibrium) absorbing state phase transition when 
the microscopic connection probability p crosses a thresh- 
old value p c = 0.6447[8]. In this article, we show that, 
instead, if sites are connected by multiple channels, the 
percolation can occur dis continuously by tuning relevant 
parameters. 

There has been a long debate on the possibility of 
discontinuous absorbing state transition in one dimen- 
sional systems. It is conjectured by Hinrichsen[9( that ab- 
sorbing state transition can not occur discontinuously in 
one dimensional systems under generic circumstances as 
strong fluctuations usually destabilize the ordered phase. 
This hypothesis challenged the existence of the discontin- 
uous transition in one dimensional triplet creation model 
observed by Dickman and Tome 1 10], and further sup- 
ported by Cardoso and Fontanari[ll|. Later studies show 
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FIG. 1: (Color online) Illustration of bond directed percola- 
tion in (1 + 1) dimension with N = 4 channels. Each site 
on a tilted square lattice has N — 4 valves (discs), of which 
n = 2 are open (blue) and the rest are closed(red) . Existence 
of at least one channel which has open valves at both ends 
(shaded), represents a bond (solid lines) between the corre- 
sponding sites. Activity can percolate only through a bond 
along a preferred direction (time) activating new sites (filled 
circles). These bonds are shown as thick blue lines. 



that the this transition is actually continuous [12], Il3j . 
Again, the apparent discontinuous transition in annihi- 
lation/fission process [14J, and the spreading process on 
a diffusing background|15j . turned out[9j to be transient 
phenomena, crossing over to a continuous transition after 
very long time. Genuine 1 st order absorbing state transi- 
tions are either seen in systems with special symmetries 
or in higher dimensions fi~7l. [l8j . 

In this Letter we introduce a multichannel model of 
bond directed percolation in (d+1) dimensions where ev- 
ery pair of neighbouring sites has N different channels 
joining them. Each channel contains a valve at every site 
that can either be open or closed independent of other 
valves. A bond is said to connect two sites if there exists 
at least one channel joining them, which has open valves 
at both the ends. Like the usual model of bond directed 
percolation, the spreading agent here propagates in a spe- 
cific direction starting from a random site. We show that 
the system undergoes a discontinuous percolation tran- 
sition in the large N limit, even in (1 + 1) dimension, 
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when the number of open valves n at each site crosses a 
threshold value y/N. We believe that these results con- 
vincingly establish the existence of first order absorbing 
state transition in (1 + 1) dimension. 

First let us discuss this multichannel model of directed 
percolation in (1 + 1) dimension in details. Its extension 
to higher dimension is straightforward. The sites on a 
tilted square lattice are connected to their neighbours by 
N directed channels (Fig. [T]). Each channel contains a 
valve at every lattice site so that, at every site, there are 
N valves in total. Each valve can either be open or closed 
independent of other valves. Flow through a channel is 
possible only when the channel is open, i.e., it has open 
valves at both ends. In analogy with the usual bond 
percolation process, a bond is said to be present between 
two neighbouring sites when there exists at least one open 
channel connecting them. The bond is absent if none of 
the connecting channels are open. At each site, randomly 
chosen n valves out of N are opened, which mimics a 
porous medium with randomly distributed directed bonds 
on (1+1) dimensional tilted square lattice. Like the usual 
bond directed percolation problem, percolation is said to 
occur here if there exists a directed path which spans the 
system. 

DP is often regarded as a dynamical spreading process 
where the fluid percolates in a given direction, interpreted 
as time. A similar interpretation for this multichannel 
model follows. In (1+1) dimension, the lattice sites in 
horizontal and vertical directions are denoted by spa- 
tial and time coordinates i = 1, 2 . . . L and t respectively 
(Fig. [1]). From any arbitrary site (i, t) the spreading fluid 
can independently propagate to the nearest neighbouring 
sites (i ± l,t + 1) only if the target site is connected to 
(i,t) by a bond. In other words, the fluid can reach (ac- 
tivate) (i ± l) th site at time t + 1 from (i, t) if there is at 
least one open channel between (i ± 1, t + 1) and (i, i). If 
none of the sites at t + 1 gets activated, i.e. if the num- 
ber of active sites at time t + 1 is zero, which happens 
when there is no bond connecting (i, t) and (i ± 1, t + 1) 
for any i, the spreading agent can not propagate further. 
Corresponding configuration at t + 1 is called absorbing. 
In this dynamical description, the system is said to be 
percolating when the fluid keeps spreading for arbitrar- 
ily long time. Note that when the system is percolating, 
there must exist one or more directed paths spanning the 
(1 + 1) dimensional lattice. Thus, this dynamic descrip- 
tion is completely equivalent with the geometric interpre- 
tation of the model discussed in the previous paragraph. 
Henceforth we continue with the dynamical description. 

A continuous phase transition from the absorbing to 
active state occurs in the usual models of bond di- 
rected percolation when the connection probability p is 
increased beyond a critical threshold p c . The steady state 
density of active sites p s plays the role of the order param- 
eter as it takes non-zero value only in the active phase 
p > p c . In the multichannel model introduced here n 
plays the role of the control parameter as the connec- 
tivity between neighbouring sites increases when more 
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FIG. 2: (Color online) (a) The steady state density of active 
sites P of usual bond directed percolation on a tilted square 
lattice with L — 1000 is plotted against p, along with 7 versus 
p [calculated from Eq. ([2])] for the mufti channel model in (1 + 
1) D with N = 1024. For any value of 7, which corresponds to 
a unique p (thus p ) one can read out p s using Q [following 
the arrow]. The resulting p s for N — 1024 and 10 12 are 
shown as solid lines in (b). Symbols therein correspond to 
p a obtained from direct numerical simulation of our model in 
(1 + 1) D with L = 1000, N = 1024. 



valves are open. Clearly there can not be any open chan- 
nel for n = 0, hence there are no bonds available for 
spreading. Corresponding steady state density of active 
sites p s = 0. Again, for n > N/2 every pair of neigh- 
bouring sites has at least one open channel connecting 
them. Starting from any random site the spreading agent 
surely activates all the sites in the steady state, resulting 
in p s = 1. Our aim here is to find if the system percolates 
for any non-zero n c < N/2. 

Undirected bond percolation on similar multichannel 
models t 19| are known to exhibit percolation transition in 
one and higher dimensions when the number of channels 
N — » 00. The transition occurs discontinuously as the 
number of open valves n crosses n c = \/_V '. In the present 
work we investigate the possibility of an absorbing state 
transition by varying n along with N as n = N 1 where 
< 7 < 1. 

Since at every lattice site, n valves are chosen randomly 
from the total N and opened, which can be done in C„ 
possible ways, the probability that there is no common 
open valve between any two neighbouring sites, i.e. the 
probability that neighbouring sites are not connected by 
any open channel, is q — C N ~ n /C™ ■ Using Stirling's 
approximation, q can be expressed in terms of v = n/N 
and N, 



q = 
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Note that, here p — 1 — q is the probability that neigh- 
bouring sites are connected by a bond. In the large N 
limit, v = TV 7-1 vanishes and g{y) — v 2 to the leading 
order in v. This results in 



p = 1 - exp(-_V 27 - 1 ). 



(2) 
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Thus, 



lim p 




for 7 < 1/2 
for 7 = 1/2 
for 7 > 1/2. 



(3) 



Clearly, the connection probability p shows a discontinu- 
ity at 7 C = | in the iV — ^ oo limit yielding a transition 
from the absorbing state (where there are no bonds) to an 
active phase (where all neighbouring sites are connected 
by bonds). Thus for any e = (7 — j c ) > the spread- 
ing agent eventually activates all the sites in steady state 
resulting in p s = 1. This absorbing state transition is 
discontinuous as the order parameter p s changes discon- 
tinuously from to 1 as 7 crosses r ) c = \- 

The order parameter p s (7;iV) for any finite N can be 
calculated from the knowledge of the steady state density 
of active sites p s DP (p) of usual bond directed percolation 
model for a given connection probability p, 



Ps (r,N) = p? p (p = i-N 2 '>- 1 ). 



(4) 



As the analytical form of pf p (p) is not known, we first 
obtain the same numerically by simulating the usual 
bond directed percolation problem on a tilted square lat- 
tice with L — 1000 for < p < 1. p s {l\ N) can be read 
out from this data using Eq. (j4j where p = 1 — N 2 " 1 ^ 1 . 
is used. This procedure is illustrated in Fig. [2ja) for 
N = 1024. The resulting p s {j;N) is shown in Fig. ^b) 
as a solid line. The symbols therein correspond to the 
same obtained from direct simulation of the multi chan- 
nel system on a tilted square lattice with L — 1000 hav- 
ing iV = 1024 channels. It is clear from this figure that 
the transition occurs discontinuously near -f c = 1/2. The 
order parameter |0 s (7;iV) for N = 10 12 , obtained us- 
ing the above procedure, is also shown in Fig. [2jb) to 
illustrate that, as expected, the transition point shifts 
towards 7 C = 1/2 as TV— > 00 . 

The existence of a discontinuity in the order param- 
eter p s is clear from the above discussions. However 
the I s * order phase transitions possess several remark- 
able dynamical aspects including hysteresis which distin- 
guish them from the continuous ones. In the following we 
will show the numerical evidence of hysteresis, the most 
important signature of a discontinuous transition. Be- 
fore going into the details of these studies first we check 
that the time-series of the instantaneous density of active 
sites p(t) shows a sharp drop when the tuning parame- 
ter 7 is suddenly decreased below the critical threshold 
7 C [marked as t — in Fig. Ela)]. Such a finite drop does 
not occur in the models of usual directed percolation as 
the order parameter there vanishes continuously at p c . 

In any discontinuous transition, it is expected that the 
system near the critical threshold fluctuates between the 
ordered and the disordered phases leading to interest- 
ing dynamical behaviour. However it is difficult to study 
these signatures in a discontinuous absorbing state tran- 
sition, as once a system reaches an absorbing configura- 
tion it remains there forever. To overcome this difficulty, 
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FIG. 3: (Color online)p(f) near the transition point: The 
system in steady state with 7 = -y c + ei is quenched to 7 = 
7c — £2 at t = and the 7 is changed back to its initial 
value after 12000 MCS. (a) and (b) correspond to our model 
with 5 = 0, 10 respectively, (c) shows p(t) for conventional 
DP near the critical point p c with same ei, ei and <5. Here 
ei = 0.005 and e 2 = 0.001. 



Bidaux et. al. [17] have introduced a special technique 
where inactive sites are spontaneously activated with a 
small rate 5. This rate S does not change the discontin- 
uous nature of the transition but converts the absorb- 
ing state into a fluctuating low density state allowing 
one to study hysteresis and related dynamical aspects 
of the discontinuous absorbing state transition. In pres- 
ence of S, the transition in this system occurs from a 
high density to a low density phase. Figure [3jb) shows 
the time series of p(t), when the system in steady state 
at 7 = 7c + ei = 0.505 (ei = 0.005) is quenched to 
7 = 7c- £2 = 0.499 (e 2 = 0.001) and kept there for 12000 
Monte Carlo steps (MCS). In the Monte-Carlo simulation 
we have used S — 10~ 5 . Here, p(t) drops to the low den- 
sity state soon after the quench (t = 0) and again jumps 
back to the high density phase as 7 reverts back to the 
initial value. This behaviour is typical of 1 st order phase 
transitions. A similar variation of the tuning parameter 
across the critical point of ordinary directed percolation 
does not show [Fig. [3Jc)] any visible change in p(t). This 
is because of the fact that in absence of the absorbing 
state, the steady state density does not change apprecia- 
bly across p c . 

Hysteresis: To study hysteresis the system is first re- 
laxed at n > VN such that the corresponding 7 is well 
above the critical threshold. Then n is decreased in steps 
of 1 after every tjt Monte-Carlo steps (MCS) until a suf- 
ficiently low value of 7 below 7 C is reached. Now the 
path is reversed i.e. n is increased in steps of 1, and 
the hysteresis cycle is completed. In Fig. [4] we have 
shown the hysteresis curves in the 7-p plot for a system 
of L = 1000, N = 1024 with different t R = 10, 30. The 
fact that the hysteresis curve bounds a finite area rein- 
forces the claim that this absorbing state transition is 
discontinuous. As expected this area decreases with the 
increase of tR. 
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FIG. 4: (Color online) The multichannel model shows hys- 
teresis near the transition point. Here activity is created 
spontaneously with a small rate 8 = 10 -5 . In the Monte- 
Carlo simulation we have used L — 1000, N = 1024. The two 
loops correspond to the relaxation times tn = 10 (circle) and 
tR = 30 (plus). Arrows indicate the direction in which the 
loops are generated. 



Site DP: The model can be recast into a different form 
by associating a set S of n integers, randomly chosen 
from a larger set {1, 2, . . . N}, to each lattice site. In this 
picture, a bond is said to connect two neighbouring sites 
(i,t) and («±l,i+l) if S(i,t)nS(i±l,t + l) ^ 0, where 
is the null set (an empty set without any element). 
The larger set, here, corresponds to N valves existing 
at each site and the smaller set S corresponds to the n 
open valves. In fact, this picture is useful in describing 
site directed percolation in this multichannel model. 

In case of site directed percolation, each site i at time t 
contains a set S(i, t) which is either a null set (inactive) 
or a set of n integers randomly chosen from {1,2,... N}. 
A site is said to be active when the associated set is 
non-empty (S ^ 0). From any given configuration at 
t the next configuration at time t + 1 can be generated 
as follows. If both S(i ± l,t) = 0, then S(i,t + 1) is 
assigned the null set 0. Otherwise, i.e. when at least 
one of them is not 0, a trial set T of n integers is created 
and S(i, t + 1) is updated as, 



5(i,t + l) 



T if T n K == 
if TOK = ®. 



(5) 



Evidently, through this dynamics, new active sites are 
created only from the existing neighbouring active sites 
with probability p which follow Eq. @ when n varies as 
n = N 1 . Thus, in N — > oo limit, here too one expects a 
discontinuous transition similar to that of the bond per- 
colation model. Details of this study are omitted here. 

It is straight forward to study this multichannel model 
of directed percolation on any kind of lattice in any (d + 
1) dimension. In higher dimensions also, one expects 
the transition to occur discontinuously at 7 C = | (i.e. 



where K is one of S(i ± 1, t) if both are non-empty; oth- 
erwise K is the only non-empty set among S(i ± 



n c = VN) as the connection probability p [Eq. ([2])] itself 
becomes discontinuous when number of channels N — > 
oo. 

In summary, we have studied directed percolation in 
(d+l) dimension on a multi channel system where neigh- 
bouring sites are connected by N channels, each having 
an operating valve at either ends. Of these A~ valves 
at each site, n are open. The fluid can percolate in a 
given direction, only through the channels which have 
open valves at both ends. In the N — > oo limit, this 
system shows a discontinuous transition from a percolat- 
ing phase to a non-percolating (absorbing) state as n is 
decreased below a threshold n c — Va". 

It is rather surprising that the absorbing state transi- 
tion in this (1 + 1) dimensional model is discontinuous. 
The possibility of such a transition has been rebutted 
by Hinrichsen Q, who argued that unless there is an 
unusually robust mechanism to prevent fluctuation from 
destabilizing the ordered phase, it is impossible to have a 
discontinuous absorbing state transition in generic (1 + 1) 
dimensional systems with short range interactions. How- 
ever, this argument does not seem to apply here, as the 
discontinuity in the order parameter is the artefact of 
the discontinuity that appears in the resulting connec- 
tion probability while varying the control parameter 7 
continuously across j c = 1/2; the role of fluctuation is 
immaterial here. 

It is worth mentioning that depinning or synchro- 
nization, which are described by a multiplicative noise 
Langevin equation[20| similar to that of DP, show dis- 
continuous transitions in certain model systems [2l|. It 
would be interesting to investigate whether the discontin- 
uous absorbing state transition obtained in the present 
study falls in the realm of their theory. 
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